AD-A145  716  MOMENT  INDUCED  BV  LIQUID  PflVLOflD  DURING  .SPIN-UP  UITHOUT  1/ 
R  CRITICAL  LAVER.  .  (U)  RRMV  ARMAMENT  RESEARCH  AND 
DEVELOPMENT  CENTER  ABERDEEN  PROVIN  .  C  H  MURPHV  AUG  84 
UNCLASSIFIED  RRBRL-TR-82581  SBI-AD-F2B8  461  F/G  19/1  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS  - 1963 -A 


AD- A 145  716 


AD/ 


TECHNICAL  REPORT  ARBRL-TR-02581 


MOMENT  INDUCED  BY  LIQUID  PAYLOAD 
DURING  SPIN-UP  WITHOUT  A 
CRITICAL  LAYER 


Charles  H.  Murphy 


August  1984 


IIS  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTEA 

BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND.  MARYLAND 


Approved  for  public  release;  distribution  unlimited. 


DT1C 

ELECTE' 
SEP  1  3  1984 

E 


S4  09  05  023 


•  «  *  *  *  v'  O  A.'  •  '  •  *  \ 


.  ^ 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Additional  copies 
from  the  National 
U.  S.  Department 
22161. 


of  this  report  may  be  obtained 
Technical  Information  Service, 
of  Commerce,  Springfield,  Virginia 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position,  unless  so  designated  by  other 
authorized  documents. 


The  use  of  trade  names  or  manufacturers'  names  in  this  report 
does  not  constitute  indorsement  of  any  commercial  product. 


Mwmlalai 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Dele  Entered) 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 


2.  GOVT  ACCESSION  NO. 


ECHNICAL  REPORT  ARBRL-TR-02581  A]) Alfa  'll!* 


4.  TITLE  (and  Subtitle) 

MOMENT  INDUCED  BY  LIQUID  PAYLOAD  DURING 
SPIN-UP  WITHOUT  A  CRITICAL  LAYER 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  6  PERIOD  COVERED 


Final 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORO) 


B.  CONTRACT  OR  GRANT  NUMBERfa) 


CHARLES  H.  MURPHY 


9-  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Ballistic  Research  Laboratory 
ATTN:  DRXBR-LFD 

Aberdeen  Proving  Ground,  Maryland  21005-5066 


ft.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

US  Army  Ballistic  Research  Laboratory 
ATTN:  DRXBR-OD-ST 

Aberdeen  Proving  Ground*  Maryland  21005-5068 


<4.  MONITORING  AGENCY  NAME  A  ADDRESS/lf  different  (rom  Controlling  Ottlco) 


16.  DISTRIBUTION  STATEMENT  (ol  title  Report) 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  6  WORK  UNIT  NUMBERS 


RDT&E  1 L 16 1 102AH43 


12.  REPORT  DATE 

August  1984 


<3.  NUMBER  OF  PAGES 

55 


IS.  SECURITY  CLASS,  (ol  t hie  report) 


UNCLASSIFIED 


15*.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  the  ebetrect  entered  in  Block  20,  It  different  from  Report) 


IB.  SUPPLEMENTARY  NOTES 


This  report  supersedes  IMR  774,  dated  April  1983. 


19.  KEY  WORDS  (Continue  on  reverea  elde  It  neceemmry  and  Identify  by  block  number) 


Critical  Layer 
Liquid-Filled  Shell 
Spin-Up 
Stabi 1 ity 


20.  ABSTRACT  rctmtbeue  an  tmrereo  WR  If  n meeeeerr  md  Identity  Or  Slock  number) 


(bja) 


For  a  fully  spun-up  liquid  payload,  eigenfrequencies  and  liquid  side 
moments  induced  by  projectile  coning  motion  have  been  computed  by  a  linear 
boundary  layer  theory.  For  partially  spun-up  liquids,  this  theory  can  fail  due 
to  the  presence  of  a  critical  layer  in  which  the  local  angular  velocity  of  the 
spinning  motion  is  near  the  angular  velocity  of  the  projectile's  coning  motion. 

(continued) 


00  /jSTti  1473  EDITION  OF  I  NOV  U  IS  OBSOLETE 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fBIVwn  Dele  Entered ) 


Ik1 


I 


SECURITY  CLASSIFICATION  OF  This  RAOEflFhan  Data  lklm« 


Late  in  the  spin-up  process,  when  more  than  90%  of  the  spin  angular  momentum 
has  been  acquired  by  the  liquid,  a  critical  layer  usually  does  not  exist.  In 
this  report,  a  linear-boundary-layer  theory  is  developed  to  predict  eigenfre- 
quencies  and  the  liquid  moment  during  late  spin-up  times.  The  predicted  eigen- 
frequencies  agree  well  with  those  computed  by  the  linearized  Navier-Stokes 
technique  of  Gerber  and  Sedney.  The  side  moments  predicted  by  the  linear¬ 
boundary-layer  theory  are  available  to  projectile  designers  for  flight 
stability  analysis. 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGEfWhan  Data  Enl.r.rf) 


TABLE  OF  CONTENTS 


Page 


LIST  OF  ILLUSTRATIONS .  5 

I.  INTRODUCTION . 7 

II.  BASIC  EQUATIONS .  9 

III.  EIGENVALUES .  17 

IV.  STEADY-STATE  SOLUTION .  18 

V.  LIQUID  MOMENT .  23 

VI.  NUMERICAL  CALCULATIONS .  24 

Table  1.  Comparison  of  Navier-Stokes  and  Linear-Boundary- 
Layer  Eigenvalues  for  m  =  0,  k  =  2,  c/a  =  .995, 

<  =  .500 .  25 

Table  2.  Comparison  of  Navier-Stokes  and  Linear-Boundary- 
Layer  Eigenvalues  for  m  =  1,  k  =  1,  c/a  =  .600, 
k  =  .443 .  25 

/II.  SUMMARY .  26 

ACKNOWLEDGMENT .  27 

REFERENCES .  37 

APPENDIX  A .  39 

LIST  OF  SYMBOLS .  43 

DISTRIBUTION  LIST . 49 


Access 

-ntis 
dtic  t 

Unanno 

Justif 

Ion  For 

GRA&I  tt 

AB  1J 

unced  □ 

By - - - - -  j 

Distribution/ _ _ 

Availability  Codes _ 

Dist 

Avail  and/or 
Special 

1 

1 


3 


LIST  OF  ILLUSTRATIONS 


Figure  Page  , 

1  Laminar  Flow  Solutions  of  Eq.  (2.5)  for  Re  =  40,000,  ! 

c/a  =  4,  <j>  =  800  rad .  28 


2  Laminar  Flow  Solutions  of  Eq.  (2.5)  for  Re  =  40,000, 

c/a  =  4,  <p  =  2600  rad. .  29 

3  Liquid  Angular  Momentum  Ratio  vs  $// Re  for  c/a  =  4 .  30 

A 

4  Spin  Angles  for  which  w  =  0.11  vs  Reynolds  number, 

c/a  =  4 . ° .  31 

5  Eigenfrequency  vs  Spin  Angle  $  for  Re  =  40,000, 

c/a  =  3.12,  k  *  3,  n  *  1,  <  -  0.443.  The  Curves  are 

from  Sedney-Gerber  NS  Theory  (Fig.  5,  Ref.  16),  the 

Circles  are  from  Linear-Boundary-Layer  Theory .  32 

6  Eigenfrequency  -r3n  vs  Spin  Angle  $  for  Re  *  5000, 

c/a  =  3.297,  k  =  3,  <  =  0.443.  The  Curves  are  from 
Sedney-Gerber  NS  Theory  (Fig.  6,  Ref.  16),  the  Circles 


are  from  Linear-Boundary-Layer  Theory .  33 

7  Liquid  Side  Moment  Coefficient  CLSM  vs  Spin  Angle  for 

j  t  =  .08  and  the  Conditions  of  Figure  5... .  34 

8  Ratio  of  Maximum  Side  Moment  Coefficients  vs  Roll  Angle, 

<t>  for  Re  =  4  x  10\  2  x  106,  c/a  =  3.12,  k  =  3,  n  =  1 .  35 

9  Ratio  of  Maximum  Side  Moment  Coefficients  vs  Roll  Angle, 

|  for  Re  =  4974,  c/a  =  3.297,  k  =  3,  n  =  1,2 .  36 

i 


j 


I.  INTRODUCTION 


The  prediction  of  the  complete  moment  exerted  by  a  spinning  liquid  pay- 
load  on  a  spinning  and  coning  projectile  has  been  a  problem  of  considerable 
interest  to  the  Army  for  some  time.  For  a  fully  spinning  liquid,  the  linear 
side  moment  was  first  computed  by  Stewartson1  for  an  inviscid  payload  by  use 
of  eigenfrequencies  determined  by  the  fineness  ratio  of  the  cylindrical  con¬ 
tainer.  Wedemeyer2  introduced  boundary  layers  on  the  walls  of  the  container 
and  was  able  to  determine  viscous  corrections  for  Stewartson's  eigenfrequen¬ 
cies,  which  could  then  be  used  in  Stewartson's  side  moment  calculation. 
Murphy3  then  completed  the  linear  boundary  layer  theory  by  including  all 
pressure  and  wall  shear  contributions  to  the  liquid-induced  side  moment.  The 
Stewartson-Wedemeyer  eigenvalue  calculations  have  been  improved  for  low 
Reynolds  numbers  by  Kitchens  et  al.4  through  the  replacement  of  the  cylindri¬ 
cal  wall  boundary  approximation  by  a  linearized  Navier-Stokes  approach.  Next 
Gerber  et  al.5”8  extended  this  linearized  NS  technique  to  compute  better  side 
moment  coefficients  for  Reynolds  numbers  less  than  10,000.  Finally,  the  roll 
moment  for  a  fully  spun-up  liquid  was  computed  by  Murphy.7 

Since  liquid  payloads  can  require  as  much  as  half  the  flight  time  to 
reach  full  spin-up,8  considerable  theoretical  effort  has  been  expended  on 

1.  K.  Stewartson,  "On  the  Stability  of  a  Spinning  Top  Containing  Liquid ," 
Journal  of  Fluid  Mechanics ,  Vol.  5,  Part  4 ,  September  1959 ,  pp.  577-592. 

2.  E.H.  Wedemeyer,  "Viscous  Correction  to  Stewartson' s  Stability  Criterion ," 
Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground,  Maryland ,  BRL 
Report  1325,  June  1966.  (AD  489687) 

3.  C.H.  Murphy,  "Angular  Motion  of  a  Spinning  Projectile  With  a  Viscous 

Liquid  Payload ,"  Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground, 
Maryland,  BRL  Memorandum  Report  ARBRL-MR-03194 ,  August  1982.  (AD 
A118676)  (See  also  Journal  of  Guidance ,  Control,  and  Dynamics,  Vol. 
6 ,  July-August  1983 ,  pp.  280-286.) 

4.  C.W.  Kitchens,  Jr.,  N.  Gerber,  and  R.  Sedney,  "Oscillations  of  a  Liquid 

in  a  Rotating  Cylinder:  Solid  Body  Rotation,"  Ballistic  Research  Labora¬ 

tory  ,  Aberdeen  Proving  Ground,  Maryland,  BRL  Technical  Report  BRL-TR- 
02081 ,  June  1978.  (AD  A057759) 

5.  N.  Gerber,  R.  Sedney,  and  J.M.  Bartos,  "Pressure  Moment  on  a  Liquid- 

Filled  Projectile:  Solid  Body  Rotation,"  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  Maryland,  BRL  Technical  Report  ARBRL-TR-02422 , 
October  1982.  (AD  A120567 ) 

6.  N.  Gerber  and  R.  Sedney,  "Moment  on  a  Liquid-Filled  Spinning  and  Nutating 
Projectile:  Solid  Body  Rotation ,"  Ballistic  Research  Laboratory ,  Aberdeen 
Proving  Ground,  Maryland,  BRL  Technical  Report  ARBRL-TR-02470 ,  February 
1983.  (AD  A125332) 

7.  C.H.  Murphy,  "Liquid  Payload  Roll  Moment  Induced  by  a  Spinning  and  Coning 
Projectile,"  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland,  BRL  Technical  Report  ARBRL-TR-02521 ,  September  1983.  (AD 
A133681)  (See  also  AIAA  Paper  83-2142,  August  1983.) 

8.  A.  Mark,  "Measurements  of  Angular  Momentum  Transfer  in  Liquid-Filled 
Projectiles ,"  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland,  BRL  Technical  Report  BRL-TR-02029 ,  November  1977.  (AD  A051056) 


predicting  the  spin-up  process  and  its  effect  on  the  liquid-induced  moment. 
Wedemeyer9  developed  a  very  simple  model  of  the  spin-up  process,  which  was 
extended  by  other  authors.10”12  The  calculation  of  eigenvalues  and  side 
moment  during  spin-up  is,  however,  a  difficult  task.  Karpov13  made  use  of  a 
Wedemeyer  suggestion  to  obtain  an  approximation  for  the  liquid  side  moment. 
Reddi 14  attempted  to  solve  the  inviscid  perturbation  equations  but  had 
considerable  difficulty  with  a  singular  critical  layer  which  is  usually 
present.  A  critique  of  this  work  is  given  by  L^nn15  who  did  get  a  solution 
for  axisymmetric  waves  for  which  no  critical  layer  is  present.  Sedney  and 
Gerber16  extended  their  linearized  Navier-Stokes  perturbation  to  calculate 
spin-up  eigenfrequencies  in  the  presence  of  a  critical  layer. 

Although  Sedney  and  Gerber  have  been  able  to  overcome  the  mathematical 
difficulties  associated  with  the  singular  critical  layer,  their  method 
requires  the  solution  of  six  complicated  first  order  differential  equations. 
This  critical  layer  occurs  in  an  annular  region  in  which  the  angular  velocity 
of  the  liquid  spinning  motion  is  near  the  angular  frequency  of  the  projec¬ 
tile's  coning  motion.  Fortunately,  late  in  the  spin-up  process  the  critical 
layer  does  not  usually  exist,  and  the  simple  linear  boundary  layer  theory  can 
be  successfully  extended  to  predict  late  spin-up  eigenfrequencies  and  liquid 
side  moment.  As  we  shall  see,  this  approach  involves  the  solution  of  only  one 
second  order  differential  equation. 


9.  E.H.  Wedemeyer ,  " The  Unsteady  Flow  Within  a  Spinning  Cylinder ,"  Journal 
of  Fluid  Mechanics,  Vol.  20,  Part  3,  1964,  pp.  383-399.  (See  also  BRL 
Report  1225,  October  1963,  (AD  43 1846).) 

10.  C.W.  Kitchens,  Jr.,  "Ekman  Compatibility  Conditions  in  Wedemeyer  Spin-Up 
Model,"  Physics  of  Fluids,  Vol.  23,  Part  S,  May  1980,  pp.  1062-1064. 

11.  C.W.  Kitchens,  Jr.,  N.  Gerber ,  and  R.  Sedney,  " Spin  Decay  of  Liquid- 

Filled  Projectiles ,"  Journal  of  Spacecraft  ana  Rockets,  Vol.  15, 
November-December  1978,  pp.  348-354.  ( See  also  BRL  Report  1996,  July 

1977 ,  (AD  A043275),  and  BRL  Report  2026,  October  1977,  (AD  AD50311).) 

12.  R.  Sedney  and  N.  Gerber,  " Viscous  Effects  in  the  Wedemeyer  Model  of  Spin- 
Up  From  Rest,"  Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground, 
Maryland,  BRL  Technical  Report  ARBRL-TR-02493 ,  June  1983.  (AD  A129506) 

13.  B.G.  Karpov,  " Dynamics  of  Liquid-Filled  Shell:  Instability  During  Spin- 
Up,"  Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground,  Maryland, 
BRL  Memorandum  Report  BRL-MR-1629 ,  January  1965.  (AD  463926) 

14.  M.M.  Reddi,  "On  the  Eigenvalues  of  Couette  Flow  in  a  Fully-Filled  Cylin¬ 
drical  Container,"  Franklin  Institute  Research  Laboratory,  Philadelphia, 
PA,  Report  No.  F-B2294,  1967. 

15.  Y.M.  Lynn,  "Free  Oscillations  of  a  Liquid  During  Spin-Up,"  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  BRL  Report  No. 
1663,  August  1973.  (AD  769710) 

16.  R.  Sedney  and  N.  Gerber ,j  "Oscillations  of  a  Liquid  in  a  Rotating  Cylinder 
Part  II.  Spin-Up ,"  Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground 
Maryland,  BRL  Technical  Report  ARBRL-TR-02489 ,  May  1983 .  (AD  A129094) 
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In  this  report  we  will  derive  this  theory  for  a  liquid  near  full  spin-up 
and  determine  validity  bounds  in  terms  of  flight  time.  Eigenfrequencies  com¬ 
puted  by  this  basically  Invlscid  theory  are  shown  to  compare  very  well  with 
those  computed  by  the  Gerber-Sedney  linearized  Navier-Stokes  theory.  Finally, 
side  moment  calculations  are  made  to  show  the  significant  changes  that  occur 
during  this  late  spin-up  period. 


II.  BASIC  EQUATIONS 


We  will  consider  a  fully-filled  cylindrical  cavity  with  radius,  a,  and 
height,  2c.  According  to  Wedemeyer,9  the  spin-up  flow  is  driven  by  Ekman 
layer  flows  at  the  flat  end  walls.  These  Ekman  layers  draw  fluid  in  near  the 
axis,  impart  rotation  to  it,  and  eject  it  near  the  cylindrical  wall.  The  flow 
exterior  to  these  Ekman  layers  on  the  end  walls  and  a  boundary  layer  ( Stewart - 
son  layer)12  on  the  cylindrical  lateral  walls  can  be  approximated  by  a  two- 
dimensional  rotating  flow  with  circumferential  velocity  VQ.  A  projectile 

attains  full  spin  during  10-20  milli-seconds  of  motion  down  a  gun  barrel. 
If  the  spin  is  assumed  to  grow  from  an  impulsive  start  at  zero  spin  to  a 

constant  value  of  I,  the  resulting  circumferential  velocity  has  the  form 


V.  -  r  4,  w  (r,  4) 


(2.1) 


a 

where  4  *  $t.  The  boundary  conditions  for  w  are* 


w  (r,0)  -  0 


0  <  r  <  a 


(2.2) 


w  (a,  4)  ■  1 


<t>  >  0 


(2.3) 


<f>  >  0 


(2.4) 


Wedemeyer  derived  a  partial  differential  equation  for  w  (r,  4)  which  was 
numerically  integrated  by  Sedney  and  Gerber. 12  The  equation  has  the  following 
form: 


■S  +  V  \zi.r»]  =  ♦  (3/r)  (M)' 

34  _  3r  Re  dr  \  3 r/ _ 


(2.5) 


*  The  proper  numerical  way  to  treat  these  discontinuous  conditions  is  given  in 
Reference  17. 


17.  R.  Sedney  and  N.  Gerber ,  " Treatment  of  the  Discontinuity  in  the  Spin-Up 
Problem  with  Impulsive  Start,”  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  Maryland,  BRL  Technical  Report  ARBRL-TR-02520 ,  September 
1983.  (AD  A133682) 


•  •  «■  f 


where 


V  =  -k£  (1  -  w) 


laminar  Ekman  layers 


=  -kt  (1  -  w)8/5  (r/a)3/5 
k£  =  k  (a/c)  (Re)"1/2; 


kt  =  .035  (a/c)  Re' 


turbulent  Ekman  layers 


k  is  usually  0.5 


Wedemeyer  derived  a  very  simple  laminar  Ekman  layer  solution  which  has  a  dis¬ 
continuity  in  shear  at  an  interior  point. 


w  =  0 


1  -  (rf/r Y 
1  -  (rf/a)' 


r  <  r 


rf  <  r  <  a 


(2.6) 


where 


rf  =  a  e 


According  to  this  solution,  the  flow  field  is  divided  into  two  regions:  an 
inner  core  flow  which  is  stationary  and  an  outer  rotating  flow.  The  boundary 
r  =  Tf  is  an  inward  moving  front  where  the  velocity  has  discontinuous  deriva¬ 
tives.  A  numerical  solution  to  Eq.  (2.5),  of  course,  does  not  have  this  dis¬ 
continuity.  rf,  however,  is  a  useful  measure  of  the  extent  of  the  spin-up 

process.  Laminar  Ekman  layer  solutions  of  Eq.  (2.5)  have  been  computed  for 
Re  =  40,000,  c/a  =  4.0  and  are  compared  for  an  early  time,  <|>  =  800  rad  (r^  = 

.6)  in  Figure  1  and  for  a  late  time,  $  =  2600  rad  (r^  =  .2)  in  Figure  2.  The 

extent  of  spin-up  can  be  described  by  the  ratio  of  the  liquid  angular  momentum 
to  the  fully  spun-up  liquid  angular  momentum: 


3  * 

r  w  dr 


(2.7) 


In  Figure  3,  I  is  plotted  versus  $/ /Re  for  several  conditions.  Note  that  for 
larger  Reynolds  numbers,  spin-up  requires  many  more  revolutions  -  i.e.,  much 


more  flight  time.  For  $  =  1000  s"1,  the  liquid  angular  momentum  is  95%  of 
that  of  a  rigid  body  in  1.6  seconds  for  laminar  Ekman  layers  (Re  =  40,000 

and  in  9  seconds  for  turbulent  Ekman  layers  (Re  =  10^). 

We  now  assume  that  velocity  components  and  pressure  can  be  expressed  in  a 
more  general  form  than  that  used  in  Reference  3. 


Vx  =  R  j  u$  e 


s 4>  -  imo 


(2.8) 


Vp  =  R  j  v$  es*  "  im0j  (a$) 


(2.9) 


V0  -  r*  w%  R{ws  es*  *  im0J 


(2.10) 


D  2 .2  n  .  o  „  4>  -  ime>  ,  a2*2v 

P  -  pL  a  $  P  +  R  psey  |(p|_a<t>) 


(2.11) 


1 

where  P  =  a"2  [w  (rj)]2  Tj  drj 

.  n 


(2.12) 


m  =  0,  ±1,  ±2  .  .  . 


*  # 

The  basic  flow,  w  (r,<|)),  is  a  function  of  time  (d>  =  it);  but  we  will  assume 
slow  variation  with  time  in  comparison  with  tne  frequency  of  the  coning 
motion. 

For  positive  spin,  the  linearized  Navier-Stokes  equations  become 
(s  -  raw  i)  v$  -  2ww$  +  a  =  Re"1  ^7Q2  v$  -  (v$  -  2  m  i  w$)  a2/r2l  (2*13) 


(s  -  m 


w  i)  w$  + 


2w  +  r  vc  -  im  a  p./r 


1  [  702  ws 


(w$  +  2  i  m  vs)  a2/r2 


(2.14) 


(2.15) 


(s  -  in  w  1)  u$  +  a^-=  Re’  vQ  u$ 


3  (r  v$) 
3r 


3us 

m  w  i  +  r  — —  =  0 
s  3x 


(2.16) 


where 


2  2 

3  .  T  m 

-  ~T  "  “2 

r3r  3x  r 


(2.17) 


Next  we  assume  the  perturbation  variables,  ug,  v$,  w$,  ps  are  not  affected  by 

the  Ekman  layers  and  Stewartson  layer  and  must  satisfy  the  boundary  conditions 
independently  by  means  of  their  own  boundary  layers.  This  assumption  is 
clearly  somewhat  arbitrary  but  allows  us  to  get  a  first  estimate  of  the  effect 
of  spin-up  flow  on  the  coning  motion.  This  will  be  refined  in  the  future  as 
our  understanding  of  the  fluid  mechanics  improves. 

The  payload  is  assumed  to  be  in  coning  motion  described  by: 

0  +  i  a  *  K  eS<f)  (2.18) 


where 


K  = 


K10e 


T  <t> 


10. 


s  =  (e  +  i)  x 


Since  (x,  r,  e)  are  cylindrical  coordinates  in  an  inertia  non-coning  axis  sys¬ 
tem,  it  is  convenient  to  introduce  cylindrical  coordinates  (x,  r,  e)  in  an 
aeroball istic  axis  system  that  cones  with  the  projectile.  The  linear  rela¬ 
tions  between  these  coordinates  are: 


x  =  x  -  R  Jr  K  es*  “  ie  J  (2.19) 


r  =  r  +  R  j  x  K  es*  '  ie  j  (2.20) 


8  =  9  -  R  ji  (x/r)  K  es*  "  19  j  (2.21) 
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In  the  aerobal 1 istic  coordinates  the  walls  of  a  cylindrical  container  whose 
center  is  at  the  projectile's  center  of  mass  have  the  simple  coordinates 

r  =  a  and  x  =  t  c.  At  these  surfaces  the  liquid  velocity  components  must 
equal  the  wall  velocity  components: 


u$  =  (s  -  i)  (r/a)  K 


vs  =  -(s  -  i)  (x/a)  K 


(m  =  1) 


(2.22) 


ws  .  [l  ♦  is  tr^Inl]  (x/a)  K 


us  "  vs  =  ws  *  0 


(m  *  1) 


(2.23) 


The  perturbation  variables  are  separated  into  inviscid  and  viscous  compo¬ 
nents  where  the  viscous  components  are  solved  by  boundary  layer  approxima¬ 
tions: 


us  ■  usi  *  usv 


V  =  v  •  +  V 
s  si  sv 


ws  =  wsi  +  wsv 


(2.24) 

(2.25) 

(2.26) 
(2.27) 


Ps  =  Psi  +  Psv  (2.27) 

For  the  fully  spun-up  liquid,  the  boundary  layer  analysis  leads  to  the  follow¬ 
ing  boundary  condition  for  tne  inviscid  flow. 


3  v  s' 
5  — si 
a  3r 


+  C 


usi 
3X  J 


=  (i-s)  (x/a)  K 

(m  *  1} 

i 

=  0 

(m  *  1) 

=  (s  -  i)  (r/a)  K 

(m  *  1) 

±c 

=  0 

(m  *  1) 

(2.28) 


(2.29) 
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(2.31) 


6C  =  -[2(m+is)J  [(2-m-is)  a  -  (2+m+is)  0  ] 

a  =  (c/a)  [s  -  (2+nOi]1/2  Re1/2  (2.32) 

3  =  (c/a)  [s  +  (2-m)i]1/2  Re1/2  (2.33) 

and  the  complex  roots  are  selected  to  have  positive  real  parts.  and  6C  are 

complex  parameters  whose  real  parts  act  like  displacement  thicknesses,  as  can 
be  seen  from  Eqs.  (2.28  -  2.29).  Their  imaginary  parts  indicate  time  lags  in 
the  fluid  response  to  the  periodic  wall  motion. 

During  spin-up,  the  lateral  wall  boundary  layer  has  the  same  Reynolds 
number  and  is  essentially  unaffected,  but  the  end  wall  boundary  layer  has  an 

exterior  flow  of  rw$  rather  than  r$.  As  a  crude  approximation,  we  will  assume 
Eq.  (2.28)  still  applies  during  spin-up  and  the  Reynolds  number  in  Eq.  (2.29) 

A 

should  be  modified  by  an  average  w. 


usi  1  (w'a>'l/2  c  sc 


3  U 


SI 


ax 


=  (s  -  i)  (r/a)  K  m  =  1 


x=±c 


=  0 


m  *  1 


(2.34) 


where 


f  • 


r  dr 


(2.35) 


The  inviscid  velocity  components  satisfy  Eqs.  (2.13  -  2.15)  for  infinite 
Reynolds  number. 


(2.37) 


■[” 


(a/r)  m  w  p$i  -  (s  -  m  w  i)  a  — ^ 


“si  *  <a/r>  ■  (s  -  m  “  "  Psi  +  (2  “  +  r  f?)  a  — 5F1]  0-1 


(2.38) 


where 


2,  ,a2  a  „  3  (*T) 


=  s  ■  2  m  w  si  +  (4  -  m  )  w  +  r 


(2.39) 


Eqs.  (2.36  -  2.38)  can  now  be  substituted  in  Eq.  (2.16)  to  yield  a  differ¬ 


ential  equation  for  ps-. 


*  r-1  fl  -  r  f  O-1]^ 

9r2  L  Sr  J  ar 

-  r-2  jjri2  +  2m  r  i  (s  -  mw  i)'1  D  9  Psi 


(2.40) 


(s  -  m  w  1)"*  D 


Next  we  make  the  usual  assumption  that  pg^  is  a  sum  of  products  of 
functions  of  r  and  functions  of  x: 


Psi  "  L  Rk  (r>  Xk  (x) 


(2.41) 


xk" ♦  c‘2  \2  xk  • 0 


(2.42) 


R.  "  +  r 


D  -  ''•l?o'1]Rk' 


(2.43) 


-  r"2  [m2  -  r2c-2  ^  +  2  m  r  i  (s  -  m  w  i)-1D  R  k  *  0 


.VLVVVj.S'.v 


where 


(2.44) 


A  O  A  O  O 

-  -  (s  -  m  w  i)  4  D  Ak 

and  Ak  is  a  constant.  If  the  constant  is  zero,  k  will  be  assigned  the  value 
of  zero,  i.e.,  aq  =  0.  The  solutions  to  Eq.  (2.42)  are 

Xo  »  A0  +  Bq  (x/c)  (2.45) 

Xk  =  Ak  cos  (a^  x/c)  +  Bk  sin  (Ak  x/c)  k  *  0,  (2.46) 

For  constant  w  and  k  *  0,  Eq.  (2.43)  becomes  Bessel's  equation.  For  a  fully- 
filled  container,  psi-  must  be  bounded  for  r  =  0  and  so  only  Bessel  functions 
of  the  first  kind  need  be  considered. 

Rk  =  CJm  (  \  r/c)  k  t  0  (2.47) 


III.  EIGENVALUES 


The  solution  for  the  perturbation  functions  can  be  written  as  the  sum  of 
a  steady-state  response  to  the  coning  motion  and  a  transient  solution  for 

which  K  =  0  in  the  boundary  conditions  (Eqs.(2.28,  2.34)).  The  steady-state 
solution  only  involves  perturbation  functions  for  m  =  1  while  the  transient 
solutions  can  have  any  integer  value  of  m  in  the  homogeneous  boundary  condi¬ 
tions.  The  homogeneous  version  of  Eq.  (2.34)  yields  a  very  simple  condition 
for  ^  when  ^  6C  is  small.  (See  page  97  of  Reference  3  for  exact  relations.) 

\  »  k  (n/2)  [1  +  (wa)“1/2  6C]  (  3.1) 


The  homogeneous  Eq.  (2.28),  however,  with  Eq.  (2.37)  imposes  a  more  compli¬ 
cated  requirement  on  Rk  (r): 


Rkv  (a,s)  -  a  Sa 


(a,s)  =  0 


(3.2) 


where* 

Rkv  (r,s)  *  [(s  -  m  w  i)  a  Rk ' (r)  -2m  (a/r)  w  i  Rk(r)]  D_1(r,s)  (3.3) 


Values  of  s  that  satisfy  Eq.  (3.2)  are  the  eigenvalues,  skn,  and  the  cor¬ 
responding  Xk  Rk  products  are  the  eigenvectors.  For  a  stable  steady-state 
solution,  these  functions  should  decay  to  zero,  which  requires  R  skn  <  0. 

A  fairly  straight-forward  iteration  can  be  used  to  find  the  skn  except 
near  singularities  of  Eq.  (2.43),  i.e.,  when 

s  =  m  w  i  (3.4) 

Since  skn  is  primarily  imaginary,  i  Tkn,  this  condition  for  the  critical 
radius,  r  =  rc,  is  approximately 

xkn  =  m  w  (r£,  *)  (3.5) 


*The  subscript  "v"  identifies  the  function  of  R^  that  is  used  in  computing 
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According  to  Eq.  (3.5),  the  critical  layer  is  located  where  the  liquid  angular 
velocity  is  t/m  and  an  internal  resonance  occurs.  For  an  inviscid  calcu¬ 
lation,  very  large  perturbation  functions  result;  and  the  much  more  complete 
linearized  Navier-Stokes  method  is  necessary.  For  any  fixed  time,  the  minimum 

A  A  A 

value  of  w  is  its  value  at  r  =  0,  i.e.,  wQ  =  w(0,4>).  For  a  large  portion  of 
the  spin-up  period,  wQ  is  zero  and  then  it  grows  to  unity,  which  indicates  a 
fully  spun-up  condition.  In  Figure  4,  spin  angles  for  which  wQ  is  0.11  are 

plotted  versus  Reynolds  number.  Since  projectile  coning  frequencies  are 
usually  less  than  0.10,  this  figure  allows  us  to  make  estimates  for  the  region 
of  validity  of  the  theory  of  this  paper.  At  a  Reynolds  number  of  40,000,  for 
example,  critical  layers  should  not  exist  at  projectile  coning  frequencies  for 
spin  angles  greater  than  1200  rad.  For  a  spin  rate  of  1000  rad/s,  this  means 
that  the  theory  of  this  report  should  apply  after  1.2  seconds  of  flight. 


IV.  STEADY-STATE  SOLUTION 

Since  the  boundary  conditions  (Eq.  (2.22))  imposed  by  the  forcing  func¬ 
tion  are  specified  for  m  =  1,  the  steady-state  solution  will  also  have  m  =  1. 
In  order  to  satisfy  these  inhomogeneous  boundary  conditions,  we  will  express 
psi  by  a  Ak  =  0  term  plus  two  series  of  products  of  functions  of  x  and  of  r. 


Psi  =  Po  +  P*+  Pe 


=  -(c/a) 


N*  Ne 


ReoXo  +  L  d*k  R*k  Xtk  +  £  den  Ren  Xen  <*> 

k  =  1  n  =  1 


(4.1) 

A 

K 


where  XQ  =  AQ  +  BQ  (x/c) 

Reo  is  the  solution  of  Eq.  (2.43)  for  A k  =  0, 


Reo  (°)  =  0;  Reo  (a)  3  1. 


eo 


X*k  =  A*k  cos  (xoV  x/c)  +  sin  U.lr  x/c)» 


'tM. 


tk 


£k 


xen  a  Aen  cos  (*enx/c>  +  Ben  sin  (xen  x/c) 


Rtk*  Ren  are  so1utions  of  Ecl*  (2*43)  for  m  =  1,  \k  =  a^ 
or  Agn,  respectively,  and  Rk ( 0)  =  0,  a  Rk * ( 0)  =  1/a. 


The  velocity  components  us^,  vs-j,  and  ws^  can  be  computed  from  Eq.  (4.1)  by 
use  of  Eqs.  (2.36  -  2.38). 


1  £ 

usi  (s  -  w  1 )_1  B0  Reo  +  c  £  d£k  Rtk  x;k 

k=l 


c  T*  d  R  X'  K 
en  en  en 


vsi  ~  Reov  Xo  +  S  d£k  R*kv  Xji 


+  V  d  R  X  K 
en  env  en 


wsi  "  (C//a)  Reow  Xo  +  E  djik  RAkw  Xi 


+  y  d  R  x  k 

en  enw  en 


where,  if  G  (r)  denotes  Reo(r),  R  k(r)  or  Ren(r),  we  have 


G  (r)  =  {  [s  -  i  w(r)]  a  G'  (r)  -  2  i  (a/r)  w  (r)  G  (r)  }  D"A(r)  • 


G„(r)  =  -  f  2w  +  r 


(4.2) 


(4.3) 


(4.4) 


G'(r)  +  i  (a/r)  (s  -  iw)  G(r)  D-1(r)  . 


The  first  series  will  be  constructed  to  satisfy  the  inhomogeneous  bound¬ 
ary  condition  on  the  lateral  wall  (Eq.  (2.28))  and  the  second  series,  the 
inhomogeneous  boundary  conditions  (Eq.  (2.34))  on  the  end  walls.  Each  series 
will  also  be  constrained  to  have  no  effect  on  the  boundary  conditions  on  the 
other  walls.  This  means  the  first  series  must  satisfy  the  homogeneous  con¬ 
ditions  at  the  end  walls  and  the  second  series  the  homogeneous  conditions  at 
the  lateral  wall. 


•*•  X^k  (±C)  +  C  (wa)  6C  X,,JLk  (*c)  =  0 


3  Renv  (a‘s) 

Renv  ^a»s)  a  6a  3?  0 


Eq.  (4.5)  imposes  simple  conditions  on  Afck  and  the  A^'s. 


A,,  =  ( nk/ 2 )  [1  +  (wj-1/2  6  ] 


A*k"° 
BAk  =  0* 


k  odd 


k  even 


Next,  we  expand  (x/c)  in  a  finite  series  in  sin  (x  .  x/c)  : 


x/c  ■  £  V  sin  <\k  x/c> 

k  =  1 


k  odd 


The  coefficients  can  be  determined  by  a  least  squares  fit.3 


£  bjk  a«k  *  bj 
k  =  1 


k,  j  odd 


where 


bj  =  c“  V  x  sin  ( A^j  x/c)  dx 


bjk  =  C"  \  Sin  x/c)  s1n  ^XAk  x/c)  dx 


X0  can  now  be  expressed  in  terms  of  sin  (A),  x/c). 


Xo(x>  ■  Ao  *  8o  £  aik  s1n  <xk  x/c) 
k  =  1 


(4.12) 


Equations  (4.1 

3  V 

% 

pute  vsl-  and 

*  “  -  * 

upsa 

cos  ( X^  x/c)  and  > 

£v 

expressions  can  be 
tion  (Eq.  (2.28)). 

;•£ 

on  A0,  and  B^ 

o 

n 

o 

< 

Equations  (4.1,  4.12)  can  be  used  in  conjunction  with  Eq.  (2.37)  to  com- 
3  v  • 

v..,-  and  — 5= —  at  r  =  a  as  a  sum  of  terms  involving  1,  sin  (xJk  x/c). 


dik  R*kv  "  a6a  3r 


A*k  =  0 


=  atkF;  B£k  *  1# 


r  =  a 


k  odd 


k  even 


(4.13) 

(4.14) 

(4.15) 


where 


F  -  1  - 


B_  R. 


The  expressions  for  p0  and  p^  are  now  considerably  simpler. 
p0  -  -B0  (x/a)  Reo  (r) 


(4.16) 


p£  =  -(c/a)  £  dik  Rik  (r)  sin  (x^  x/c)  k  =  1,3,5...  (4.17) 

k  =  1 


The  construction  of  the  pe  series  and  the  determination  of  B0  is  somewhat 
more  complicated.  The  primary  difficulty  is  the  determination  of  the  Xen's. 
These  are  the  eigenvalues  that  satisfy  Eq.  (4.3).  The  subscript  n  identifies 


the  eigenvalues  in  order  of  increasing  size.  These  eigenvalues  can  be  found 
by  an  iterative  process.  The  first  guess  for  these  can  be  obtained  from  a 
table  of  inviscid  eigenfrequencies  such  as  Table  3  of  Reference  3.  In  this 
table,  the  eigenfrequencies  t  ,  are  tabulated  as  functions  of  a  reduced 
fineness  ratio,  f*. 


Tnk  ‘  Fn  <f*> 


(4.18) 


where  f*  =  c/ka  =  irc/2x^ 


(4.19) 


These  tables  are  constructed  by  finding  pairs  of  t^'s  and  x^'s  that  satisfy 
Eq.  (4.3)  for  Sa  =  0.  Thus,  for  a  specific  frequency  t  =  Tn((,  each  of  the  Fn 
functions  can  specify  a  wave  number  Xgn  =  X^  that  satisfies  Eq.  (2.43).  For 
example,  for  t  =  .04,  and  c/a  =  3,  xel  =  4.524,  xg2  *  9.364,  xg3  *  14.405. 

Equation  (4.2)  can  now  be  used  in  boundary  conditions  (2.34)  to  yield 


Aon  =  0 
en 


B  =  1 
en 


n  *  0 


(4.20) 


1  R0„  +  c  y  d  \  [cos  x  +  (wj-1/^  5  x  „  sin  x  1  R„„ 
o  eo  n  -  1  en  en  en  'a'  a  en  en  en 


=  (s  -  i)2  g(r) 


where  g(r)  =  (r/a) 


(4.21) 


Next  we  fit  (g(r)  to  a  series  in  the  Ren‘s  by  use  of  least  squares. 


9(0  =  E  R 


en  en 


n  =  0 


(4.22) 


where 


Lc .  a  =  c . 
jn  en  j 


n  =  0 


j  =  0,1,. ..N, 


(4.23) 


c .  =  a 


9  ( r )  Ki  dr 


(4.24) 


(4.25) 


S' 


Rej  Ren 


dr 


Eqs.  (4.21-4.22)  can  now  be  combined  to  give  conditions  on  Bn  and  den. 


B0  =  (s  -  i)2  a 


eo 


(4.26) 


dar,  t°os  Kn  +  (wj-1/2  <5  A  „  sin  A  1  =  (s  -  i2  a  n  *  0  (4.27) 

en  en  en  a  c  en  en  en 


For  a  fully  spun-up  liquid,  Re0  is  r/a,  aeo  =  1,  aen  =  0  (n  *  0),  pe  =  0,  and 

the  solution  reduces  to  that  of  Reference  3.  When  the  liquid  is  not  fully 
spun-up,  the  aen's  are  not  zero  but  are  moderately  small.  Only  when  the 

expression  in  brackets  is  small  does  dgn  become  large.  This  is  the  case 
when  xgn  is  near  a  x^. 

V.  LIQUID  MOMENT 

The  linear  liquid  moment  in  response  to  the  coning  motion  described  by 
Eq.  (2.13)  can  be  expressed  in  the  following  form  in  terms  of  the  non-spinning 

aeroballistic  axis  system,  X,Y,Z  which  pitches  and  yaws  with  the  missile. 

MLY  +  1  \z  "  V2*2-1  (CLSM  +  1  CLIm)  *  e$*  (5.1) 

The  imaginary  part  of  the  dimensionless  moment  coefficient,  CLIM,  causes  a 
rotation  in  the  plane  of  the  coning  moment  and  is  called  the  liquid  in-plane 
moment  coefficient  while  the  real  part,  Cl$m,  causes  a  rotation  out  of  this 
plane  and  is  called  the  liquid  side  moment  coefficient.  Since  Cl$m  affects 
the  damping  of  the  coning  motion,  and  only  affects  the  frequency,  pre¬ 
diction  of  Cl5M  is  the  primary  objective  of  any  theory. 


The  major  components  of  the  liquid  moment  are  due  to  the  pressure  on  the 
lateral  and  end  walls  of  the  container.  Lesser  components  are  due  to  the 
viscous  wall  shear  on  the  lateral  and  end  walls;  thus,  the  liquid  moment 
coefficient  can  be  given  as  the  sum  of  four  terms: 

T  <CISM  *  1  CIIN>  *  V  *  "pe  *  "vt  *  "ve  (5'2) 

The  pressure  as  given  by  Eq.  (2.11)  Is  specified  as  a  function  of  r  and 
x.  Eqs.  (2.19  -  2.20^  can  be  used  to  yield  a  linear  approximation  for  pres¬ 
sure  as  a  function  of  r  and  x. 


23 


%  %  %  •.  %\\ 


p  (pL|a2)-1  -  p  (r)  +  ^g-  (r  -  r)  +  R  j  p$es*  ‘  10J 


»S  /V  *2 

v  r»  A" 


Mr)  *R  j[Ps 


(5.3) 


The  pressure  moments  relative  to  the  center  of  the  cylinder  are: 


.  C  „  2  TT 


.  .  m 


-  i(2nacK)_1  e_s^  j  x  e1 9  R  j  [ps  -  (x/a)  K  j  es< 


Sc 

x  fps  (a,x)  K"1  -  (i/a)  ]  dx 


dx  de 


(5.4) 


Sa 

[ps  (r,c)  K"1  -  rw2  (c/a2)]r2  dr  (5.5) 

n 


The  viscous  components  of  the  liquid  moment  are  computed  by  the  formula 
(7.2,  7.3)  of  Reference  3  and  are  given  in  Appendix  A. 

VI.  NUMERICAL  CALCULATIONS 

In  Reference  16,  Navier-Stokes  calculations  of  eigenfrequencies  are  given 
for  m  =  0,  Re  =  43,000  and  for  m  =  1,  Re  *  5000,  30,000,  40,000,  50,000, 

2  x  10®.  6C  was  taken  to  be  zero,  and  laminar  Ekman  layers  were  assumed  for 

all  cases  except  for  Re  of  2  *  10®  which  involved  turbulent  Ekman  layers.  The 
rotational  symmetric  case  with  m  *  0  is  given  in  Table  1.  The  degree  of  spin- 
up  Is  indicated  by  r^/a  as  defined  in  Eq.  (2.6).  This  parameter  varies  from  1 

to  0  as  time  varies  from  0  to  »  (fully  spun-up).  For  this  case,  there  is  no 
critical  layer;  and  the  simple  Linear-Boundary-Layer  theory  (LBL)  gives  the 
same  frequencies  as  the  numerically  more  difficult  Navier-Stokes  theory  (NS). 
The  damping  rates  do  differ  by  .0005.  Since  the  damping  rates  are  determined 
by  the  viscosity  on  the  cylindrical  wall,  this  indicates  the  inaccuracy  of  the 
boundary  layer  approximation  at  this  relatively  low  Reynolds  number. 

Table  2  gives  results  for  the  rotational ly  asymmetric  case  of  m  =  1.  In 
the  first  part  of  this  table,  the  eigenfrequencies  for  n  »  1  are  negative  and, 
hence,  no  critical  layer  exists.  Once  again,  the  frequencies  agree  well  al¬ 
though  the  comparison  of  the  damping  rates  is  much  poorer.  In  the  second  part 

of  the  table,  for  n  »  2,  a  critical  layer  does  exist  for  r^/a  =  .27  (w  =  .02) 


r 


but  the  nonzero  value  of  seems  to  be  sufficient  to  allow  a  good  evaluation 
of  r.  The  very  small  value  of  e  does  indicate  the  nearness  of  an  instability. 
The  NS  value  of  damping  in  the  presence  of  a  critical  layer  is  a  healthy  value 
and  shows  the  power  of  this  method  in  the  presence  of  a  critical  layer. 

In  Figure  5,  the  ei genf requencies  are  compared  for  both  laminar  and 
turbulent  Ekman  layers  and  agreement  is  good  for  later  spin-up.  The  Navier- 
Stokes  values  are,  of  course,  the  only  ones  appropriate  for  early  spin-up. 

Finally,  Figure  6  shows  another  comparison  of  eigenfrequencies  for  the  very 
low  Reynolds  number  of  4974  and  the  later  spin-up  values  are  in  good  agree¬ 
ment. 

The  L8L  theory  is  the  first  to  compute  moment  coefficients  during  spin- 
up.  This  is  not  an  intrinsic  limitation  of  the  NS  theory  and  approximate 
numerical  codes  are  presently  being  developed  for  the  Navier-Stokes  theory. 
An  example  of  the  variation  of  the  side  moment  coefficient  with  time  is  shown 
in  Figure  7  for  the  conditions  of  the  laminar  Ekman  layers  curve  of  Figure 
5.  For  the  fixed  coning  frequency  of  .08,  the  side  moment  coefficient  is 

computed  for  4>  ranging  from  1000  to  2000  rad.  This  covers  the  region  between 
the  disappearance  of  the  critical  layer  and  the  establishment  of  a  steady- 
state  spin  profile.  As  could  have  been  predicted  from  Figure  5,  the  peak 
occurs  at  <j>  =  1400  rad.  This  however,  is  the  first  calculation  of  the  value 
of  this  peak  side  moment  coefficient. 

For  the  condition  of  Figure  7,  the  fully  spun-up  values  of  and  its 

associated  maximum  CLSM  are  .040  and  .312.  According  to  Figures  5  and  7,  at 

4>  =  1400,  these  quantities  are  approximately  .08  and  .36.  Thus,  the  peak  side 
moment  coefficient  at  this  point  of  spin-up  is  16%  greater  than  its  fully 

spun-up  value. 

In  Figure  8,  we  plot  the  ratio  of  this  maximum  side  moment  coefficient  to 
its  value  for  <j>  =  ®  versus  <{>.  For  the  smallest  value  of  <j>  for  which  the 
critical  layer  is  not  present,  i.e.,  1,000,  the  maximum  side  moment  coeffi¬ 
cient  is  40  %  larger  than  its  fully  spun-up  value.  In  this  figure,  the 

corresponding  plot  for  Re  =  2  x  106  is  also  given.  For  the  higher  Reynolds 
number,  the  spin-up  time  is  much  larger  and  is  in  excess  of  40,000.  At  the 
disappearance  of  the  critical  layer  for  41  of  7,000,  the  peak  side  moment 

coefficient  is  30%  less  than  its  fully  spun-up  value.  Finally  in  Figure  9, 

similar  plots  of  the  maximum  side  moment  coefficient  ratios  are  given  for  the 
conditions  of  Figure  6. 

VII.  SUMMARY 

1.  During  the  later  portion  of  the  spin-up  process,  the  period  of  the 

circumferential  flow  is  a  function  of  radial  distance  from  the  spin  axis  but 

is  never  near  the  period  of  the  coning  motion,  i.e.,  no  critical  layer  exists. 

2.  For  this  time  period,  the  simple  linear  boundary  layer  theory  can  be 
extended  to  compute  eigenfrequencies  and  side  moment  coefficients. 


3.  The  eigen frequencies  have  been  compared  with  those  computed  by  the 
much  more  complicated  Navier-Stokes  perturbation  method  and  excellent  agree¬ 
ment  has  been  obtained. 

4.  Side  moment  coefficients  have  been  computed  for  critical  layer  free 
portion  of  the  spin-up  process  and  the  maximum  values  corresponding  to  a  par¬ 
ticular  eigenfrequency  determined. 
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Figure  1.  Laminar  Flow  Solutions  of  Eq.  (2.5)  for  Re 
c/a  =  4,  4>  =  800  rad. 
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ximum  Side  Moment  Coefficients  vs  Roll  Angle 
>00,  c/a  =  3.297,  k  =  3,  n  =  1.2. 
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APPENDIX  A 


The  viscous  components  of  the  perturbation  variables  and  liquid  moment 
during  spin-up  differ  from  the  expressions  given  in  Reference  7  in  two  ways: 
(1)  the  boundary  condition  on  ws  as  given  by  Eq.  (2.22),  and  (2)  the 

replacement  of  Re  by  wg  Re  in  the  end-wall  boundary  layer  profiles. 


More  speci f ical ly ,  the  boundary  layer  profiles  on  the  lateral  wall 
(r  =  a)  are: 


w 


sv 


1  +  is  +  a 


3  w  (a) 1 
ar  J 


(x/a)  K  -  wg. 


(r-a)/a« 

a 


(Al) 


(r-a)/afi 

usv  =  -  C(i  -  s)  K  +  u$.  (a,x)]  e 


(A2) 


The  viscous  pressure  perturbation  at  the  lateral  wall  is 


Psv  (a.x)  =  2Sa  Kl  +  is  +  a  (x/a)  K  -  w s-  (a,x)  |  (A3) 


On  the  end  walls  (x  =  ±c),  a  single  complex  velocity  profile  is  required  for 
the  viscous  liquid  moment  integral: 


wsv  -  1  vsv  =  *|wsi  c)  -  1  vsi  (F»  c) 

(A4) 

-  ^2  (1  +  is)  +  r  j  (c/a)  K  je  8  c 


3  =  i  (c/a)  wa1/2  6a_1  yj  (1  -  i s ) / ( 1  +  is)  (A5) 


The  viscous  pressure  perturbation  is  zero  on  the  end  walls. 


LIST  OF  SYMBOLS 


Aq  coefficient  in  the  expansion  of  X0;  see  Eqs.  (2.45)  and 

(4.12-13) 

Ak*  Aen»  A^k  coefficients  in  the  expansions  of  Xen,  X^ 

a  radius  of  a  liquid-filled  cylindrical  cavity 

a  coefficients  in  the  expansion  (4.18)  of  r/a;  the  solution  of 

system  (4.19) 

a  .  coefficients  in  the  expansion  (4.8)  of  x/c;  the  solution  of 

**  system  (4.9) 

BQ  coefficient  in  the  expansion  of  XQ;  see  Eqs.  (2.45)  and  (4.25) 

Bk*  Ben*  Btk  coefficients  in  the  expansions  of  Xk,  Xea,  X^ 
bj  defined  by  Eq.  (4.10) 

bjk  defined  by  Eq.  (4.11) 

CLjm  liquid  in-plane  moment  coefficient 

CLSM  liquid  side  moment  coefficient 

c  one-half  the  length  of  the  liquid-filled  cylindrical  cavity 

Cj  defined  by  Eq.  (4.24) 

Cjn  defined  by  Eq.  (4.25) 

0  0  (r,t)  defined  by  Eq.  (2.39) 

den»  d*k  defined  by  Eq.  (4.1) 
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LIST  OF  SYMBOLS  (Continued) 
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m 


mpe» 


mve*  mv«. 


Ne»  na 


p0»  Pe»  P* 


Ps 

P$i»  Psv 


liquid  angular  momentum  ratio,  Eq.  (2.7) 

Bessel  function  of  the  first  kind  of  order  m 
Km  exp  (i  4>10) 

magnitude  of  the  coning  motion 
k  (a/c)  Re~*^ 

.035  (a/c)  Re"1/5 

azimuthal  wave  number  in  Eqs.  (2.8  -  2.11);  m  =  0,  ±1, 

±2,  ... 

I m  | 

non-dimensional  liquid  moment  due  to  end-wall  and  lateral 
pressure,  respectively 

non-dimensional  liquid  moment  due  to  end-wall  and  lateral 
viscous  force,  respectively 

maximum  values  of  k  and  n,  respectively,  considered  in  the 
end-wall  and  lateral  expansions,  Eq.  (4.1) 

liquid  pressure 

components  of  psi,  Eq.  (4.1) 

liquid  pressure  perturbation 

inviscid  and  viscous  components  of  p$ 


LIST  OF  SYMBOLS  (Continued) 


defined  by  Eq.  (2.12) 

radial  coordinates  in  the  inertial  and  aeroballistic 
systems,  respectively 

the  front:  the  r  value  separating  the  rotating  and  non 
rotating  fluid 

a  $/v,  Reynolds  number 

functions  of  r  in  the  expansion  of  ps1-,  Eq.  (4.1) 
defined  by  Eq.  (4.3) 
defined  by  Eq.  (4.4) 

function  of  r  in  the  expansion  for  psi,  Eq.  (2.41) 
defined  by  Eq.  (3.3) 

(e  +  i)  t 

eigenvalue  of  s:  a  value  of  s  that  satisfies  Eq.  (3.2) 
time 

components  of  the  liquid  velocity  perturbation  in  the 
inertial  system 

inviscid  part  of  us,  v$,  ws 


viscous  part  of  u  ,  v  ,  w 


LIST  OF  SYMBOLS  (Continued) 


\*  V 

y;> 


velocity  components  of  the  liquid  in  the 
inertial  cylindrical  system 

V0  (r$)-1  for  K  =  0 

b 

A 

average  value  of  w,  Eqs.  (2.34  -  2.35) 

axial  coordinates  in  the  inertial  and  aeroballistic 
systems,  respectively 

functions  of  x  in  the  expansion  of  ps^,  Eq.  (4.1) 

function  of  x  in  the  expansion  of  psj,  Eq.  (2.41) 

angles  of  attack  and  side-slip,  respectively,  in  the 
aeroballistic  system 

functions  of  (c/a.  Re,  m,  s),  defined  by  Eqs. 

(2.32  -  2.33) 

exponential  damping  per  cycle  of  coning  motion 
parameter  in  k^,  usually  0.5 

azimuthal  coordinates  in  the  inertial  and  aeroballistic 
systems,  respectively 

a  constant  in  Eqs.  (2.42)  and  (2.46) 

[-D/(s  -  m  w  1)2]l/2  Ak 

constants  in  the  expressions  for  Xgn,  ,  respectively 
kinematic  viscosity 
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LIST  OF  SYMBOLS  (Continued) 
liquid  density 

ratio  of  coning  frequency  to  spin  frequency 
eigenvalue  of  x 
|t,  spin  angle 

a 

orientation  angle  of  K 
spin  rate 
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